1. IntroductionIn quantum mechanics, the whole quantum system can be totally described when the density operator is known. Considering the abstract of the density operator, some quasi-probability distribution functions, such as Wigner function, Q function as well as P function, are introduced in phase space,[1–3] which is also an appropriate stage connecting classical and quantum mechanics. With the aid of these distribution functions, one can describe the noncalssical properties which are not involved in the classical mechanics. For example, the negative Wigner function is an obvious indicator of nonclassicality.[4] In addition, it is also convenient for expectation calculation of quantum mechanics operators.
Although the operator is abstract, the operator problem is common in quantum physics. For instance, operators’ ordering is actually an important topic in quantum optics and quantum field theory, where arranging operators are needed for path integral quantization, and converting quantum master equation (ME) to Fokker–Planck equation in quantum optics, and so on. Furthermore, it will also be beneficial for deeply understanding the essence of the operator from the view point of operator. It is an excellent example that superoperator method is introduced to deal with the dissipative problems.[5] Then these MEs are solved to obtain density operators’ infinite operator-sum representation via the thermo-entangled state representation.[6,7] In addition, the operator method has also been applied to the fields of quantum information and preparation of quantum states. For example, an operator describing the quantum-state transfer is derived such that the fidelity of quantum teleportation is thus shown to be a direct consequence of the information obtained in the measurement.[8] Recently, this operator method is used to discuss the preparation of nonclassical states via conditional measurement,[9,10] such as photon-addition/subtraction and the superposition of both, as well as quantum catalysis.
In many cases, density operator is used rather than distribution functions. For instance, quantum entanglement (negativity, entropy) can be calculated by the trace norm of partial transpose of density operator.[11] And the minimum error probability in quantum illumination can be evaluated by two possible output density operators. On one hand, there exist an evolution process for quantum state and an interaction between quantum system and its surroundings. Thus the density operator is usually difficult to obtain. On the other hand, however, it will be relatively easy to describe the resulted state by Wigner function. Thus it is a natural idea how we can obtain the density operator from the Wigner function. In this paper, we shall construct the relation between the Wigner function and density operator by introducing the Wigner operator. Then we pay our attention to the evolution of quantum states through dissipative channel.
The rest of this paper is arranged as follows. In Section 2, we review the Wigner operator and discuss its properties. Based on the property, the evolution of the Wigner function after unitary transformation is considered in Section 3, especially how to obtain density operator from the Wigner function. We use the method above to examine the evolution formula of the Wigner function in dissipative channel in Section 4. Section 5 is devoted to deriving the density operator of the TMSV in different dissipative channels. Our conclusions are drawn in the last section.
2. Wigner operator and its property2.1. Introduction of Wigner operatorFirst, we here consider the Wigner function of the density operator. For any single-mode density operator ρ, the Wigner function in coordinate representation
is defined as in Refs. [1]–[3]
Here,
W(
q,
p) is actually a quasi-probability distribution function satisfying
and its marginal distributions, the probability distribution of position and momentum, can be given by performing the integration over the variable
p or
q, respectively,
where
is the eigenvalue of the momentum operator.
In order to discuss the quantum system from the view point of operator, it is often convenient to rewrite Eq. (1) as the following form
by introducing the operator
, named as Wigner operator,
[12] whose form is given by
From Eq. (
6) it is ready to see that the Wigner operator has the following properties: (i) the Wigner operator is Hermite,
i.e.,
, which protects the real property of the Wigner function; (ii) The trace of the Wigner operator is a constant:
and for two same modes, the trace of product is just a Dirac delta function,
i.e.,
where
and the overlap relation
is used.
2.2. Normal ordering form of Wigner operatorOn the other hand, one can gain a deep understand if the specific form of the Wigner operator is obtained. For this purpose, we can appeal to the technique of integration within an ordered product (IWOP) of operators,[13–15] and use the coordinate state
whose explicit form is given by
where
and
and
a satisfying
are creation and annilihation operators, respectively. With the help of the IWOP technique, it is ready to obtain
or
where the symbol
denotes the normal ordering, and
and
are the coordinate and momentum operators, satisfying the commutative relation
. Equation (
10) is just the normal ordering form of the Wigner operator, which presents a concise Gaussian type within the normal ordering.
2.3. Weyl ordering and its invarianceActually, except the normal ordering and anti-normal ordering form of operators, the Weyl ordering form is a special interesting one. Then what is Weyl ordering form of the Wigner operator? Using Eq. (11) and the formula converting the normal ordering form to the Weyl ordering form,[16] i.e.,
where
is the coherent state, it is ready to see that the Weyl ordering form is just the following delta function, a concise form,
Here Diracʼs
δ-operator function within Weyl ordering symbol
can bring much convenience for discussing some problems in quantum optics and quantum statistics.
In addition, there is a very useful property for these operators within Weyl ordering symbol, that the Weyl ordering is invariant under the similar transformations,[17] which means
This implies that
S operator can run across the symbol
. This property can be checked by using the Weyl expression of operator or the definition of Diracʼs delta function for Eq. (
13). These properties above shall be effective for dealing with the unitary evolution systems and are going to be used in the following calculation.
3. Evolution of Wigner function after unitary transformationFrom Eq. (5) one can see that the Wigner function can be formally derived from the Wigner operator. Then how can we obtain the density operator if we know the Wigner function? For any single-mode density operator ρ, combing Eqs. (5) and (8), it is ready to have
or
Here we should notice that
. Actually,
is just the Weyl correspondence of the density operator
ρ.
[18] Equation (
15) or equation (
16) implies that the density operator
ρ can be derived from the Wigner function if the Wigner function is easily obtained. Thus equation (
15) construct the bridge between the Wigner function in phase space and quantum density operator.
Now, we consider a single mode density operator
going through a unitary evolution characteristics of unitary operator U with the following transformation,
where
,
required by unitary property, or
where
and
. Actually, for many quantum evolution, the transformation relation of operators are relatively easily obtained. Thus it is convenient to examine the evolved Wigner function from the view point of operator transformation. Under the transformation shown in Eq. (
17) and using Eq. (
15), then the evolved density operator, denoted as
ρout, can be put into the form
Using Eqs. (
13) and (
14) it is ready to gain
then substituting Eq. (
20) into Eq. (
19) yields the output density operator,
i.e.,
Comparing Eq. (
15) with Eq. (
21), it is clear to see that the Wigner function after evolution can be directly obtained by replacing
with the classical canonical transform
in the initial Wigner function in phase space,
i.e.,
Although this conclusion is drawn in the single-mode case, it is also valid for general multi-mode case. For ease of use, for multi-mode case, we show the corresponding formula as follows:
and
where
, and
. These disussions above tell us that when the matrix describing a unitary evolution shown in Eq. (
17) is known, then the evolved Wigner function can be directly obtained. According to Eq. (
15) or Eq. (
21) one can obtain the evolved density operator. That is to say, once the initial Wigner function and the unitary transformation are known, the output state and the Wigner function can be obtained by a replacement and an integration, respectively. From the above discussions, one can clearly see that it is all because of the Weyl ordering form of the Wigner operator and its main invariance property mentioned above in Eq. (
14), which is also the reason that we consider the relation between the Wigner function and the density operator. In addition, the Wigner function is a real function and can be measured experimentally.
4. Evolution formula of Wigner function in dissipative channelNext, we use the method above to examine the evolution formula of the Wigner function in dissipative channel. In most cases, in order to consider the dissipative effect including photon loss (amplitude dissipation), thermal channel and phase sensitive channel, the beam splitter models are widely used, accompanied by two inputs: signal state and accessory state, denoted as ρin and ρm respectively. Thus the output state can be expressed as
, where
and
which leads to the transformation relation
where
,
are the transmissivity and reflectivity, respectively. From Eq. (
25), one can see that
where
,
and
,
.
In order to realize our purpose by using Eq. (22), one needs to know the initial Wigner function of input states. For any two inputs ρin and ρm which can be chosen to be vacuum, thermal state, etc. here, according to Eq. (23), the corresponding Wigner function of ρout can be given by
where
. According to Eqs. (
24) and (
7), the output
is calculated as
Then using Eqs. (
5) and (
8) one can obtain the Wigner function of
Equations (
29) and (
30) can be used to derive the density operator and the Wigner function of the output when going through a dissipative channel modeled by beam splitter, respectively. Actually, equations (
29) and (
30) are general expressions for deriving density operator and Wigner function for any two inputs states when they are combined by a unitary transformation martix.
Now, we use Eq. (30) to derive the evolution of the Wigner function for two kinds of channels: photon loss and thermal channels. For photon loss channel,
, whose Wigner function is
. Substituting it into Eq. (30) one can obtain
After making a variable substitution,
, and
which leads to
, then equation (
31) can be reformed as
Equation (
32) indicates that one can derive the output Wigner function after photon loss only when the initial Wigner function is known. If we take
, equation (
32) has the same form as the result in Ref. [
19].
When the channel is a thermal one, i.e.,
where
is the average photon in the thermal state. The corresponding Wigner function is given by
In a similar way to derive Eq. (
32), by replacing
R with
, it is ready to obtain
or
Equations (
32) and (
35) are just the time evolution formula of the Wigner function when going through photon loss and thermal channels, respectively.
[19] Here we should mention that there is a Fourier transformation relation between the Wigner function and the characteristic function, thus it will be interesting to bridge the relation between the characteristic function and the transformation, which can also be further extended to multimode case.
5. Density operator of the TMSV in dissipative channelAs an often used entangled resource, the two-mode squeezed vacuum (TMSV) is widely applied in many fields, such as quantum teleportation, quantum key distribution, quantum metrology, and so on.[20–23] Here we pay our attention to the Wigner function and the density operator after the TMSV goes through dissipative channels. It will be useful for further discussion in realistic case.
Theoretically, TMSVs can be shown as
, where
is the two-mode squeezed operator with r being a squeezing parameter, and
is two-mode Fock states. Its Wigner function is
![](cpb_28_8_084206/cpb_28_8_084206_eqn37.gif) | |
5.1. TMSV through two symmetrical photon-loss channelsHere for simplification, we assume that the two-mode of TMSV go through the same photon loss channel. According to Eq. (32) the Wigner function of dissipative TMSV can be obtained as
Performing directly the integration above, one can obtain
where
It is interesting to notice that the Wigner function of output is similar to the initial Wigner function, which indicates that the Gaussianity of TMSV can be kept after photon loss. In other word, the loss channel can be seen as a Gaussian channel, although the output state is no more the pure state. In particular, when
T = 1 corresponding to the case without photon loss, equation (
39) just reduces to Eq. (
37) as expected. Further using Eq. (
15) and the normal ordering form of the Wigner operator shown in Eq. (
10) as well as the IWOP technique,
[13–15] one can directly obtain the Wigner operator after the dissipation by taking
n = 2 in Eq. (
24),
i.e.,
where
Or rewriting Eq. (
40) as
which is a Gaussian form, as expected.
[24] Equation (
41) is just the normally ordering form of the output state, which is convenient for our further discussions about the normalization, nonclassical properties including the degree of entanglement, and quantum steering.
5.2. TMSV through single-side thermal channelNext, we further examine the case that one of the TAMV goes through the single-side thermal channel. In fact, this asymmetrical case can be often encountered in quantum information. For instance, quantum illumination can be seen as this model, in which one mode of signal state is going to through a thermal channel. In the process of dealing with the problem, a beam splitter with thermal state as an input is often used. Here we take the TMSV as the inputs and one mode through the thermal channel. Under this case, using Eqs. (35) and (37), the Wigner function of output state can be derived as
where
Here we should notice that only one mode goes through the thermal, thus the parameters for other mode in phase space is kept unchanged. Similarly, substituting Eq. (
42) into Eq. (
24) we finally obtain the density operator after the thermal channel,
i.e.,
It is obvious that when
T = 1, equations (
42) and (
43) just reduce to Eq. (
37) and the normal ordering form of the TMSV, respectively. In particular, when
corresponding to the single-side photon loss,
[25] the results can be shown in Eqs. (
42) and (
43). Again, equation (
43) is the normally ordering form when one mode of the TMSV through a thermal channel. Actually, for a convenience of our calculation, we often need some different ordering forms such as Weyl ordering form, anti-normal ordering form, and
Q −
P ordering form, the normally ordering form above can be further converted into other ordering forms by using corresponding formula not shown here.
6. ConclusionsIn summary, from the coordinate representation of the Wigner function, we introduced the Wigner operator, which leads to that the Wigner function can be seen as an average value of the Wigner operator within the density operator (the considered quantum state). Then using the explicit form of the coordinate state and the IWOP technique, we derived the normal ordering form of the Wigner operator, and the Weyl ordering form. It is shown that the Wigner operator presents a Gaussian form within normal ordering, while its Weyl ordering form is just a Diracʼs delta function. These properties are very useful for discussing the evolution of quantum state from the view point of operator.
Based on these properties above, we further construct the relationship between initial state and output state after a unitary evolution, including Wigner function and density operator, by which it will be convenient and effective to obtain the output density operator. In addition, using this method, we further derived two evolution formula of the Wigner function when any single-mode quantum state goes through photon loss and thermal channels. As applications, we finally considered the two-mode squeezed vacuum going through photon loss and thermal channels, and the output Wigner function and the normal ordering form of output density operator are obtained. These results can provide some convenience for further discussion in other cases, such as quantum steering and quantum metrology as well as quantum cryptography.[26,27]
It should be emphasized that the above method can be extended to multi-mode system although we only considered single- and two-mode cases. Recently, the non-Gaussian operation and conditional measurement are used to improve the quantum correlation, such as the quantum teleportation, the degree of entanglement, quantum steering, and so on.[28] Our method is also valid for non-Gaussian states especially for obtaining the final density operator.[29] In the whole process, the IWOP technique has an important role.